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Digital Adaptive Flight 
Controller Development 

■by 

Howard Kauftnan*, G. Alag**, 
P. Berry**, and. S. Kotob** 


Summary 

Digital fly-hy-wire flight control systems are currently of interest 
■because simple mechanical linkages are no longer a"ble to cope with the many 
control problems associated with high performance aircraft and space vehicles 
In view of this, adaptive flight control algorithms have been developed for 
ultimate implementation and flight testing in an onboard digital computer. 

In order to achieve this objective, the following hierarchy of activity 
was pursued: 

. Development of a direct digital controller amenable 
to on-line adaptation. 

. A study of on-line discrete parameter identification procedures. 

. A development of digital adaptive flight control logic suitable 
for inflight implementation. 

The above analysis was performed using equations of motion and the 
related parameters for a typical fighter aircraft as supplied by NASA. 

As a result of con5)utational experiments performed using typical measure 
ment noise characteristics, two configurations are recommended, namely: 

* Assist. Professor of Systems Engineering Division, Rensselaer Polytechnic 
Institute, Troy, N. Y. 12l8l 

** Gi^iduate studerits at Rensselaer Polytechnic Institute, Troy, N. Y. 12l8l 




An interfacing of linear quadratic optimal regulator logic with 
a weighted least squares estimator. 

An interfacing of control logic, designed, using single stage 
performance indices, with a weighted least squares estimator. 



1.0 INTRODUCTION 

1.1 Background 

Digital fly-by-vd.re flight control systems are currently of interest 

because single mechanical lihkeiges are not able to cope with the many control 

12 

problems associated with hi^ perfoimance aircraft and space vehicles. ^ 

DigiteJ. implementation is extremely advantageous because of: 

The significant weight and volvime savings. 

. The availability of low cost integrated circuits. 

. The ability to design con 5 >leac controllers which previously 
were inqoossible to implement onboard sm aircraft. 

. The high reliability of digital logic. 

. The capability for time sharing multiple control loops. 

Purtheimore the need for an adaptive control system capable of control adjust- 

3 4 

ment has been established for : 

Pioviding uniform stability and handling qualities over the , 
con 5 >lete fli^t envelope despite dreistic changes in the 
open loop characteristics of the aircraft. 

Providing acceptable flying qualities over a wide range 
of external disturbances due to atmospheric turbulence and 
outer loop command signals. 

I 

Design of a digital adaptive flight control system requires careful 

consideration of the following factors; 

Computer storage limitations 
Con 5 >uter cycle time 
Computer operation time 
Stability requirements 
Parameter changes 
Control limitations 
Pilot input signals 



Control hysteresis 
Test signal requirements 
Turbulence effects 

Development of such a system can be guided by assigning a heirarchy of 
structure to the system as shown in fig. 1 . At the lowest level, directly 
interacting with the aircraft dynamics, is the Direct Digital Controller 
(DDC). This consists of a set of gains that multiply the appropriate states 
along with the necessary summing logic. Because this DDC controller operates 
directly on the aircraft, it is executed at the feistest rate. Typical san^jle 
rates might be 30 times per second. 

Because any given set of DDC gains cannot be satisfactory over this 
con^jlete flight envelope, it is necessary to incorporate some sort of logic 
that adjusts these gains whenever necessary. Thus the adaptation logic is 
placed at the second and third control levels. 

For small parameter deviations, and subsequently small deviations 
in performance, an in5>licit adaptation algorithm might be used. This does not 
require explicit parameter identification and could be used at a sample rate 
intermediate to the DDC logic and the e3q>licit adaptation logic. 

Finally, in order to adjust the DDC gains in the presence of large 
parameter variations, some type of on-line identification appears necessary 
in order to guarantee system stability. The execution rate of this logic is 
clearly dependent upon the rate of parameter variation and must be determined 
through experimentation. 

A summary of the state of the art in adaptation can be found in 
ref. 5. Some of the more salient conclusions of this study are: 

For nonminimum phase plants, the only methods that can 
achieve good control are those using some form of 
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FIG. 1 Digital adaptive flight control system. 
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explicit parameter identification. 

. Adaptive techniques requiring explicit plant identifi- 
cation tend to require the most conputational capability. 

. The best adaptive control over output performance character- 
istics is potentially provided by those adaptive methods that 
explicitly identify the plant. 

. Excessive bandvd.dth (gain) is undesirable because of the 
resulting sensitivity to measurement noise and the danger 
of exciting high order structural or sensor modes. 

. Control methods which explicitly identify the aircraft have 
the lowest gain level requirements because the gains can be 
adjusted directly to their proper values which are functions 
of the airframe parameters. 

( 

Despite the previously cited studies, an actual digital implementa- 
tion of a completely adaptive flight controller has to date not been attempted. 
This may be attributed to the limitations of the available analog equipment as 
well as the reluctance to change over to digital control. Consequently, a 
research effort was initiated by NASA Langley Research Center for the purpose 
of applying modem digital control theory towards the development of an imple- 
men table digital adaptive flight controller. 

1.2 Objectives 

In view of the attractiveness of digital adaptive flight control, 
an effort to develop such a system has been pursued by the Systems Division of 


6 



Rensselaer Polytechnic Institute (RPl) since Serptember 1972. The objectives 
of this effort have been: 

1) To investigate adaptive control procedures which might 
be capable of being implemented in a digital flight 
computer. 

,2) To test these procedures using linearized lateral and 
■ longitudinal data as supplied by NASA for a typical 
fighter aircraft operating over a given flight envelope. 

3) To validate the design using one of NASA's hybrid 
simulations. 

1.3 Scope and Outline 

Development of an adaptive control system requires consideration of 

1) The system equations of motion. 

2) The selection of an index of perfoimance. 

3) The coinparison of actual process behavior with the 
desired or Ideal behavior. 

U) The parameter gain adjustment procedure used to drive 
the process more towards the desired behavior. 

Whereas in the initial phases of the development, the linearized 
latereil and longitudinal equations of motion are being considered separately, 
for control purposes, final results wall include the performance resulting 
from the application of various developed controllers to a nonlinear six- 
degree-of-freedom simulation. A description of the linearized equations is 
contained in Section 2.1. 

Because it is convenient to express desired aircraft response to 
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pilot commands in terms of a mathematical model, all performance indices con- 
sidered in the development penalize some measure of the error between the 
process and model outputs. A detailed outline of such "model following" 
performance is foimd in Section 2.2. 

Usage of such a model also facilitates the instantaneous con 5 >arison 
of process performance with the desired behavior. This leads to the possi- 
bility of implementing either in^jlicit adaptation logic which does not require 
process parameter identification or explicit adaptation logic which indeed 
does use online parameter estimates. Following a discussion of these alter- 
natives in Section 2.3, various control structures amenable to adaptation are 
presented in Section 2.4. Procedures for online identification are then 
described in Section 2.5* 

Finally based upon the various subsystem designs given in Section 
2.0, overall adaptive control logic is presented in Section 3 » 0 , and con^ta- 
tional res'ults based on this logic are given in Section 4.0. 

Recommendations for implementation and for further studies are 
found in Section 5*0. 

1.4 Significance 

Development of a digital adaptive flight control system is of 
significance not only to the particular aircraft considered but also to digital 
process control in general. Such a development r^resents an important appli- 
cation of modem digital control theory that is a step towards narrowing of 
the gap between theory and practice. The proposed logics while feasible for 
digital process computers, is clearly impractical for analog implementation. 

Of immediate significance, however, is the demonstration that 
digital fly-by-wire technology is capable of providing desirable handling 



qualities over the con5)lete flight envelope of a high perfoiroance aircraft 
despite changes in the open loop characteristics of the airframe. 
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2.0 REQPIREMEDIJTS FOR DIGITAL ADAPTIVE FLIGHT CONTROL 

Items essential to adaptive control design are individually discussed 
in the following four sections. System r^resentation is first developed in 
Section 2.1 followed by a discussion of performance evaluation in Section 2.2. 
Adaptation principles and control structure are presented in Sections 2.3 and 
2.4 respectively. Finally^ methods for online parameter identification are 
presented in Section 2.5. 

Restilts pertinent to the tuning and the performance of these procedures 
are given in Section 4. 

2.1 System Representation 

The linearized dynamics of the aircraft as supplied by NASA may be 
r^resented by the vector state equation 

• 

X = F X + G u . . (2,1) 

where x denotes the state vector 

u denotes the control vector 

and F and G are matrices of the appropriate dimensions. 


For linearized lateral notion 



roll rate 
yaw rate 
sideslip angle 
roll angle 


and u = , 

IK] 

1 aileron deflection' 

• 1 

1 

1 rudder deflection 


/ 
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For linearized longitudinal motion 


X = 

lA 

^pitch rate \ 


A - 

velocity j 


i“ 

angle of attack I 

. ■ 

w 

^pitch angle / 


and u = ((j = (elevator deflection) 

The elements of F €md G, known to vary with mach number and 
jJ^titude, were provided for a typical filter aircraft for the six flight 
conditions (FC's) shown in the fli^t envelope of Figure 2. These are provided 
in Appendix A. 

The objective of the research was to find implementable digitaLL 
algorithms for conqputing the control signals S , ^ , euad h so as to 

£t 6 - 

insure iiniform and desirable handling quantities for an aircraft flying within 
the given flight envelope. This was to be performed assuming that during flight 
the elements of F and G were not readily available (e.g., as scheduled 
functions of mach number and altitude). 

Because of the need to inq)lement the control system in a digital 
flight con 5 >uter as shown in Fig. 3, Eqn. 2.1 was transformed into the equiva- 
lent discrete form: 

X (k+l) = A X (k) + B u (k) (2.2) 

FT 

where A = e 

(C “ 

T = Sampling Time 

and x(k), u(k) denote x, u at time k.»T. 
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EXCESSIVE SPEED SLUGGISH 
ADVERSE YAW INSTABILITY IN PITCH 



12 


FIG. 2 Flight Envelope 








FIG. 3 Digital Flight Control Loop 









Eqn. 2.2 is a valid representation asstiming that F and G do not vary for 
k T ^ t ^ (k+l) T and that the control signal is constant between sample 

times, i.e., 

u (t) = u (kT) for kT ^ t ^ (k+l) T 

This latter relation will he true assuming the use of a zero order hold type 
of digital to analog converter. The foimer assumption is valid if the rates 
of change of F and G are such that changes over the sample period are 
negligible. 

2.2 Performance Evaluation ^ 

Inherent to the effectiveness of any adaptive control system is the 
capability for rapidly assessing the performance and making the necessary 
modifications to the control gains. One such procedure that fits these 
requirements and at the same time has the potential for insuring uniform 
handling qualities is the concept of model following control as depicted in 
fig. 4. This concept has been of interest to many investigators over the 
past few years. ’ In fact, relative to these efforts, Erzberger has 
published a set of "perfect model following" conditions under which the out- 
put of the process can be made identical to the output of the model. 

Being that the ideal objective of model following flight control is 
to force the aircraft to respond as the model would to a given pilot command,' ' 
it is often desirable to simulate the online model dynamics in the flight 
computer and to generate the aircraft control signal using the actual aircraft 
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FIG. 4 Model Following Control System 










states, the pilot input cOTsnands, and the model states. This situation is 
sometimes referred to as the pilot's flying the computer while the con^iuter 
is flying the aircraft. 

More precisely, the model following problem can he stated as follows: 
Given the aircraft dynamics: 


X (k+l) =. A X (k) + B u (k) 
p' ^ P ^ P P 


(2.3) 


where; ^p(^) aircra,ft (nxl) state vector at sample time k 

Up(k) is the (mxl) control vector 


and 


A^, Bp are matrices with the appropriate dimensions ; 


find the control u (k) such that the process state vector x (k) 

P P . ■ 

approximates "reasonably well" some model's state vector Xj^^(k) defined by 


the equation: 

"m W 

where x^ (k) is the (nxl) model state vector 
Uj^ (k) is the (mxl) pilot input and 

Bjj^ are matrices of the appropriate dimensions. 

One possible variation is to find u^- (k) such that the (rxl) 

dimensional process output vector 

H X (k) 

P P 

approximates the model's (rxl) output vector 


(2.4) 


H X (k) 

m m ' ' 
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A farther possibility to consider is the situation in which there 
is proces.s measurement noise n(k) such that the only available process out- 
put is: 

=: (k) + (k) (2.5) 

Tarb-ulence can be considered in the above system by including a 
process noise term w^ such that the plant equation becomes: 

Xp (k+1) = -A^ (k) + Bp lip (k) + P\ .(2.6) 

where 

Wj^ is a stochastic process whose statistics are chosen to 
approximate the atmospheric conditions and 
P is the corresponding distribution matrix. 

2.3 A<iaptation Logic 

In the very early phases of the project, considerable study was 
given to the tradeoffs involved in selecting either an explicit or an implicit 
adaptation algorithm, the differences being that 

, In explicit adaptation^ on-line estimates of the aircraft 
parameters are used for gain adjustment. 

In implicit adaptation, some measure of the error between the 
actual and the desired state trajectories is used for gain 
adjustment. That is, no e:q)licit parameter identification is 
used. 

This study and a similar one made by TASC in 1970^ indicated that 
explicit adaptation would initially be preferable because: 
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stability was to be required in the presence of severe 
parameter variations. 

. Effectiveness of iinplementing iinplicit adaptation is 
related to the satisfaction of the conditions of perfect 
model foUowing.^^^^'^ These in general will not be 
satisfied for HASA supplied models . 

Implicit adaptive contiol con 5 )arable to the hyperstability approach 
used by Landau^ is to be considered at a later date. 

2.4 Control Algorithms 

2.4.1 Perfect Model Following 

In a model following (mf ) adaptive control system the design 
goal is to force the con^ensated system to duplicate the performance of a 
reference model. Ihe knowledge of plant parameters miables controller gains 
to be set so as to achieve the desired characteristics. 

There are two configurations of mf, known as inqolicit mf 
and real inf. As fig. 5 show^ in inplicit mf, the aircraft controls are 
formed from the aircraft states and pilot input. No dynamic coupling exists 
between the model states and the closed-loop plant; the model state x^ 
appears only in the performance index. 

Figure 6 illustrates real model following in which the model 
states must be generated for use in forming the control law. Alignment of 
plant and model in the presence of \ancertainties such as unknown parameters 
and random disturbances, requires this type of control. This enables a con- 
tinuous correction of the errors between model and plant states even in the 
presence of unknown disturbances. 
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Erzberger has shown the existence of a control for perfect 


model following provided that: 


(I - Bp b/ ) (A^ - V = 0 


( 2 . 7 ) 


(I - B B 
P P 


) B 


m 


= 0 


( 2 . 8 ) 


where B is the pseudo inverse of B . 
P P 

then the control input 


If Ihese conditions are satisfied. 


u = B ^ (a. - A. ) X + B^ B u (2.9) 

p p ' m -p p m m \ y/ 

will result in perfect model following. 

Chan^*^ has given an alternative control law which yields 

perfect following if eqs. 2.7 and 2.8 are satisfied, and which yields 

stability in the sense of boundedness otherwise. 

This controller has the form: 

u 
P 

where 

e 

controllers is the single stage design approach used by Motyka’^'” for 
simultaneously weighting both the state error and the contfol effort at a 
single time instant. Elimination of the control effort weighting and a slight 
modification of the error state yields Erzberger’s perfect mf control law. 


Ui . U2 


= K 


= B 


t 


(A-A)x + B B u 
' m p -ffl p mm 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


X - X , the error vector and K is a gain matrix chosen 
m p 

to stabilize (Ap - B^ K). 

Related to the development of perfect model following 

IP 
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2.4.2 Single Stage Optimization 

In view of the work by Chan and Motyka, and the attractive- 
ness of control gains readily adjustable in the presence of parameter changes, 
single stage indices were considered for digital adaptive control design. 

The essential feature of a single stage performance index is 
the fact that only a single point in time is considered. If e (k) describes 

the error in mf at sample time k, then a performance index of the form 

e (k) Q e (k), where Q is a positive semn definite weighting matrix, 
allows the designer to specify the relative in^jortance of mf in different 
states by trading off the size of mf errors in different states. In 
addition, the plant control signal can also be weighted in the performance 
index by including a factor of the form u^ (k) R u^ (k), where u^ is the 

control signal and R is a positive semi-definite weighting matrix. Such a 

factor allows further tradeoff between mf in the plant states and the magni- 
tude of the control signals required to effect that mf . Thus, a performance 
index of the form 


J(k) = e'^ (k) Q e(k) + u^^ (k) R u^' (k) 


(2.13) 


was considered. 


To find the vector u to minimize J, V J is formed 

— p u 

^ -p 

and set to zero. The resulting equation gives the vector u . 

R 

The first formulation, called RF for "response feedback", 

12 

as developed by Motyka for continuous time systems is now discussed here for 
discrete time. For known present plant state x (k), the next state x (k+l) 
can be required to be that state which would result if Xp(^) were used in 2.4 
instead of 2^(k) i.e., if the plant had the dynamics of the model. Such a 


21 



desired correspondence would le realized if the plant, state at (k+l), i.e., 
x^(k+l) were equal to: 

\ Xj, (k) + (k) (2.14) 

Thus the error signal e in eq. 2.13 160011168: 


(a X (k) + B u (k)) - (a X (k) + B u (k)) 

PP PP P 


p 


is set equal to zero to obtain: 
u (k) = pR + B Q B ~1 

p ' p ^ p 


-1 T 

B Q * 
P 


[~(A - a ) X (k) + B u (k)n 

•— m P P mm — • 


= K X (k) + K u (k) 

V “ 

p m 


(2.15) 


(2.16) 


The features of this contipller include the following: 

(1) It is closed-loop containing feedback only from the 
plant and not from the model. The model matrices 
are, however, used in forming the controller gains; 

(2) Steady state behavior of the plant is not necessarily 
equal to that of the model since model states are not 
used for comparison; 

(3) Real-time simulation of -the model is not necessary. 

A second fonnulation, to be denoted the "same state" algorithm 


requires that the plant states be exactly equal to the model states at the 
next sampling instant i.e., it requires that: 
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Xp (k+1) (k+1) 


(2.17) 


Substituting from (2.3) and (2.4), the error signal 

e = (A^ Xp (k) + Bp Up (k)) - (A^ x^ (k) + u^ (k)) 


(2.18) 


Again setting 'O 


= 0 gives : 


Up(k) =•[>-. Bp% Bp] - Bp 


(2.19) 


or 


u (k) = K X (k) - K X (k) + K u (k) 

p xm' X p u 

^ m P ni 


This SS controller has the foUovd.ng features: 

(1) It is closed-loop with respect to the plant and 
63cplicit with respect to the model; 

. -(2) . The factor Z = is common 

“‘v to all 'three gain matrices. 

(3) Because x^ is explicitly used, steady state 
differences between plant and model states 

might be removed given enough independent controlling 
elements . 

(4) Gains can easily be updated on-line without having to 

solve any nonlinear algebraic (Riccati) equations. 

Because of these attractive characteristics, the same state 
controller was chosen for implementation in an adaptive control loop, details 
of which are presented in Section 3.1. 
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2.4.3 Infinite Time Optimal Quadratic Control 

\ 

2. 4. 3.1 Performance 

The performance index which is selected for minimi- 
zation in 5 )lies the controller type, the degree of model- following, and the 
overall control effort expended. A typical perfoimance index might consist 
of Some positive semi-definite function of the model following error 

(x (k) - X (k)) balanced against a positive definite function of the control 
p m 

(k). Inclusion of control rate and/or integrals of the model following 

error will result in a controller with a dynamic structure, i.e., with 

7 

integral feedback terms. Asseo has demonstrated that such dynamic feedback 

can yield reduced sensitivity to plant parameter variations, and other 

13 l4 15 

investigators have examined various specific structures. 

7 

In parti cTilar, Asseo considered the control rate 

• 

Up(t) rather than the control itself as the independent variable. This 
necessitated the treatment of the plant control s-s an additional state 

with the state equation 


u (t) = V (t) 

P P 


( 2 . 20 ) 


In the discrete case, the integrator is replaced by 
an accumulator (i.e., a "unit delay with unity feedback), and the corresponding 
equation becomes : * ' - 


u^(k+l) = Up(k) + Vp(k) 


( 2 . 21 ) 


An . attractive feature of this formulation is the ability to include in the 

performance index both positive definite functions of n^(k) (to prevent 

overly large controls) and of v (k) (to prevent overly large control rates). 

P 

This control structure was the actual one utilized in the experimental studies 
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A second dynamic feedback structure as siiggested 
13 

by Tiroshi and Elliott is analogous to tbe classical proportional-integral 
feedback structure. By integrating the model-plant error, it is expected 
that any step change in the model states vdll result in zero steady-state 
plant-model error. Again the discrete equivalent of an integrator is an 
accumulator and the state equations for the "integral" are: 

X (k+1) = X (k) + C(x (k) - X (k)) (2.22) 

oi oi ^ m 

The matrix C specifies which model-following errors will be accumulated. Care 
must be taken in the selection of the matrix C, for an m controllable system 
can resTolt. A description of this difficulty and a test for system controll- 

l4 

ability is given by Porter. This controller would be especially suitable 
for a transport type aircraft (because of the accurate steady-state control) 
but may be more responsive to noisy feedback than the first structure because 
the proportional-integral controller does include direct feedback of the plant 
states themselves. 

The various possibilities suggested above were 
broken into two classes of performance indices. The performance index 
corresponding to the first controller includes terms weighting both model- 
following error and control rate, and has the form; 

QD 

1 ^ Q(Xp(k) - \W) + Vp'^(k) R Vp(k) (2.23) 

^ k=0 

where Q and R are matinees to be chosen by the designer. 

The performance index corresponding to the second 
structure includes terms weighting model-following error (^p " accumulated 
model- following error x as defined in eq. 2.22, and the magnitude of the 

, d« 
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control Up. This index^ which therefore includes three weighting matrices 
S, and R, is of the form: 

OD 

I ^ (Xp(k)-x^(k))'^ Q(Xp(k)-x^(k))+ S xjk) 

k=0 

+ Up^ (k) R u^(k) - (2i24) 

Notice that if the matrix S is all zeros, then no weight is placed on the 
acctunulated plant-model error and this controller will no longer effectively 
have any feedback dynamics. 


2. 4. 3 . 2 Controller Structure 

Infinite time optimal linear regulator theory was 
used to find the minimum of the performance indices constrained by the appro- 
priate state equations . This resTolts in a control formed by feeding back all 

15 

states through constant gains. (Details are presented in the next section.) 
With this in mind, the specific controller structure inqolied by these perfom- 
ance indices may be described. Further, the performance indices and state 
equations for both cases may be put in terms of the infinite time optimal 
linear regulator problem: 

. Find u(k) to minimize 


1 

2 


I 


k=0 


x^(k) Q’ x(k) + u"^(k) R' u(k) 


(2.25a) 


subject to 

x(k+l) = A x(k) + B u(k). 


(2.25b) 
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For the first perfoimance index which penalizes 
control rate rather than control magnitude (2.23), fig-. 7 illustrates the 
final closed-loop structure. Note that any control input (from the model 
state and control), or any plant state feedback signal must he "filtered" 


thiough the dynamic feedback controller before it effects the plant control. 
Keshaping the state equations into a form to which the optimal linear regu- 
lator results are directly applicable (2.25a and 2.25b) gives: 



Note that the” pilot's input (equal to the model control u^(k)) is modeled as 
a constant - an ansuitq>tion that does not overly distort the reality of the 
situation and allows a complete analysis of the problem from a theoretical 
viewpoint. 

The performance index that includes the accumulated 
model-following error (eq. 2. 2U) leads to a controUfer structure analogous to 
the classical proportional-integral controller, as shown in fig. 8. Manip- 
\ilating the state equations into a form to which the optimal linear regulator 
results are directly applicable, (2.25a and 2.25b) yields; 
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Xp(t) = Fp(t)Xp(t) + Gp(t) Up(t) 

Up(k + I) = Vp(k) + Up(k) 

Vp(k) = K2|Xp(k) + K22Up(k) + K23X^(k) + K24Um(k) 

FIG. 7 Dynamic Feedback System Using "Control Rate" Feedback 
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FIG .8 DYNAMIC FEEDBACK SYSTEM USING"INTEGRAL ERROR' 
FEEDBACK 


Xp(t) = Fp(t) Xp(t)+Gp(t)Up(t) 

Xq (k+i)= Xa(k) + C(Xp(k)-Xm(k)) 

Up(k)=K„ Xp(k)+K,2 Xa(k) + K,3 Xm(k) + K,4 Um(k) 


Xm(k) 
















As noted above, if the weighting matrix S were all 
zeros, then the feedback gain ^ would be zero and the controller 

would be non-dynamic . 



constrained by the state equation: 

x(k+l) = A x(k) + B u(k) (2.25b) 

15 

is, well known and derived in Moore and Anderson among others. The control 
has the constant feedback structure: , 

u(k) = - K x(k) (2.28) 

where . r 

K = (b'^PB + R)“^ b'^PA (2.29) 
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and P is specified "by the steady state Riccati equation: 

P = Q + A^PA - A^PB(R + B^PB)‘^ B^PA (2.30) 

This nonlinear equation in P may he solved in a nuniber of ways. The 
simplest appears to be by iteration of the time varying Ricatti equation 

P(k+1) = Q + A^ P(k) A - A^ P(k) B(R + B*^ P(k) B>"^ B^P(k) A (2.31) 

This equation converges to the solution for any positive definite initializa- 
tion of P(k) provided that the pair !Ia,b 3 is controllable and the pair 
[a. e] is observable (where Q = D^D). At this point it should be noted 
that even though the model is not controllable, the solution will still exist 
because of the controllability of the plant. 

The adaptation logic for this type controller is 
based upon an online iteration of eq. 2.31. Since the aircraft parameters 
vary continuously and relatively little within the anticipated gain update 
cycles, it can be expected that the exact solution to the corresponding steady 
state Riccati equation will not vary significantly between gain updates. Thus 
if at each gain update time, the Riccati equation is initialized with the most 
recent solution, it is hypothesized that it will be necessary to iterate only 
a few times to find the proper solution. 

Finally as a practical matter, it should be noted 
that' coD 5 )utation can be simplified if eq. 2 . 3 I is partitioned taking into 
account the symmetry of P and the sparseness of the aiogmented system matrices 
A and B as defined by 2.26 and 2 . 27 . The resulting expressions can be found 
in eqs. 3 . 1 , 3 . 2 , and 3 . 3 . 


31 



2.5 Parameter Identification 
2.5.1 Introduction 

Selection of a parameter identification algorithm for an 
adaptive flight control system requires consideration of the foUovdi^ items: 

. Required frequency of peirameter updating 

. Measurement noise characteristics’ 

. The desirability of estimating true stability 

derivatives versus the elements appearing in the ■ 
discrete transition matrix (These will be nonlinear 
functions of the stability derivatives . ) 

. Nonlinearities induced by simultaneous estimation 
of parameter and states 

. Performance under closed loop control 

If parameters vary relatively slowly, then a batch process 
( 17 ) 

such as quasilinearization^ might be feasible. However, if updating must 
be done more frequently, then an online or recursive procedure with some sort 
of fading memory should be employed. 

In designing parameter identification algorithms, attention 
must be given to the measurement noise characteristics. Whereas a relatively 
large variance necessitates the use of an identifier with a long memory to 
achieve smoothing, a small variance will enable the use of a short memory 
Identifier that will be more responsive to parameter variations . 

Further, it is conjectured that in a digital environment, 
greater success would be achieved if the unknown parameters of the discrete 
transition matrix itself rather than the physical stability derivatives are 
identified. This follows because the discrete transition matrix is a highly 
nonlinear function of the stability derivatives making identification rather 
difficult. Furthermore computation of the discrete control law requires 
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estimates for elements of the discrete transition matrix, and not of the 
stability derivatives. A precise description of the discrete parameter esti- 
mation problem is presented in Section 2.5-2. 

In designing an identifier, consideration must be given to 
algorithms for estimating only parameters and algorithms for estimating both 
parameters and states simultaneously. Whereas the foimer set of algorithms 
utilizes the measurements to represent the true state, the latter set depends 
upon the augmentation of the unknown parameters to the system state vector. 

As shown in ref. (l8), the colHbined state and parameter estimation problem is 
nonlinear because of the need to determine quantities that multiply each other 
(parameter times state). Consequently one must select a priming trajectory 
about which the system is to be linearized. 

lypical algorithms for parameter estimation only include the. 
weighted least squares approach and the stochastic gradient method. These are 
discussed in Section 2.5;3. Section 2.5-^ discusses two approaches for esti- 
mating both parameters and states; namely, the extended Kalman filter and a 
decoupling pi?;)cess in which the linear state and parameter estimation problems 
are treated separately and alternately. 

Finally, of importance is the performance of these identi- 
fication procedures \mder closed-loop control. Because such control often 
resTrlts in transient behavior for only a very small amount of time and steady 
state, behavior for a relatively large amount of time, there may not be 
sufficient excitation to allow accurate enou^ identification. Thus the need 
for an induced dither signal must be examined. Such dither could in fact be 
produced by feeding back. for control computation the noisy state measurements 
themselves rather than filtered state estimates. In any event the performance 
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of the identifier should ultimately he measured hy the overall behavior of 
the adaptive control system rather than the individual estimates it produces 
for the various parameters. This follows from the fact that not all states 
are equally excited by any given pilot input command. Hence various motions 
will be completely decoupled or very insensitive to the specific values 
obtained for some of the parameters. 

Thus while there may not be sufficient excitation present to 
accurately track a parameter, its value may not be very influential in deter- 
mining the control signal needeid for the maneuver being undertaken. Accurate 
tracking of al 1 parameter continuously will only be possible if dither can be 
acceptably introduced into the motion of the aircraft. 

Results of identifying under closed- loop, control using the 
weighted least squares algorithm are presented in Section 4.5. 

2 . 5.2 Problem Statement , 

The problem considered to date is that of determining the 
values of certain parameters appearing in the discretized aircraft equations 
of motions given exact measurements of the inputs and noisy measurements of the 
outputs. Noisy input measurements are to be considered in subsequent studies. 
Evaluation of the accuracy of this identification is to be based upon the 
resulting adaptive control loop performance. * 

As given by eq. 2.3, the lateral or longitudinal motion of 
the aircraft is represented by the vector difference equation 

Xp(k+ 1 ) = Ap(q) Xp(k) + (q) u^(k) ( 2 . 32 ) 

where q is used to denote a vector vdiose elements are unknown parameters 
appearing in the pleint matrices. For estimation purposes, it will be assumed 
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that the system measurements can be described by 


yp(k) = H Xp(k) + ^ (k) (2.33) 

where: (k) is a zero mean measurement noise sequence with covariance 


matrix R 


R(i,0) =(r 


where X .=1, if i=j, and S . .=0 otherwise 
XO 10 10 (2.3U) 


and: H is a selector matrix indicating o^st which states or combinations 

of states are measured. For the work considered to date, a.n states were 
assumed measurable, and hence H=I, the identity matrix. 

4 

i^rpical noise characteristics are included in Table I; these 
were provided by NASA - Langley Research Center. Although the specified noise 
is correlated, many of the experiments conducted to date have utilized 
uncorrelated (i.e., white) noise sequences, having the values cited in 
Table I, in order to sinplify the analysis. 


2.5.3 Procedures for Parameter Estimation Only 

2. 5 . 3.1 wei^ted Least Sqtiares Parameter Identification 

The simple weighted least squares linear parameter 

20 

estimation scheme was chosen to avoid some of the problems of algorithms 
which estimate both parameters and states. Since noisy states are fed back 
through the control loop and used for parameter estimation, it is expected 
that some sort of state smoother will ultimately be necessary. 

To apply this method to the systems defined by 
eqs. 2.32 and 2.33^ it is necessary to differentiate between the set of para- 


meters (in both 


and B ) that are to be identified and the set of para- 
■ P 


meters that are not to be identified. 
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NOISE CHARACTERISTICS 
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In particular a convenient representation for (2.32) 
results if all the parameters (both a (i,j) and b (i,j))which are to be 
identified are placed in a vector q , and a1 1 those parameters which are not 
to be identified are placed in a vector £ . With these definitions, the state 
equations can be rewritten as: 


x(k) = C(k. - 1) q + D(k - 1) s 


( 2 . 35 ) 


S"t 

where C^^(k-l) contain the state or control measurement at the (k-l) -sample 
that couples x.(k) to q. , and D (k-l) .contains the state or control 


measurement at the (k-l)^^ sample : that couples x.(k) to s.. A zero entry 

for a particular C. . (or D. .) would indicate that no coupling exists between 

■ 

X. and' q. (or between x. and s.). 

1. -3 ,v 1- - .a 1 

Defining the estimate fomed on the basis of .k 

samples for the parameter vector as q(k), the weighted least squares algorithm 

20 21 

as developed by Yo\ang and Mendel minimizes the performance index: ... 


J = J (y(i) - x(i))^ R"^(y(i) - x(i) ) 


( 2 . 36 ) 


1=0 


where y(i) is the actual measurement made of x(i) 

^ A 

and , , x(i) = C(i-l) q(k) + D(i-l) _s 

A ■ • ^ 

i.e., x('i) is'.^the "estimate" of x(i) formed by substituting q(k) 


into the state equation (2.35). 


-1 


The weighting matrix R is usually selected to be 


a diagonal matrix whose elements r^^ reflect either the confidence and/or 

the relative magnitudes of each of the measurements y^ . Furthermore a fading 

2i -1 

factor of the form a (a>l) can.be included in R in order to give major 
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emphasis to the most recent measurements. 

20 

Minimization of (2.36) yields 

q(k) = P(k) B(k) 

~ k n 

where P (k) = \ ( Jl-l)) R"^ C(.^5^-1) 

J 

k 

and B (k) = ^ (X-1) R’^ [y(J^) - D ( - l) s] 

1=1 

Note that computation of the estimate q. (k) from these foimula requires at 

each new sample a matrix inversion of order equal to the number of parameters 

being estimated. This solution however can be rewritten in a recursive manner 

that takes into account the effects of an additional measurement in a manner so 

as to reduce the matrix inversion requirements. In particular, the estimate q(k) , 

based upon k samples, can be written in terms of q(k-l) , the estimate based upon 

20 

k-1 samples using the equations: 

q(k) = q(k-l) - K(k) (C(k-l) q(k-l) - z(k) ) (2.37a) 

z(k) = y(k) - D(k-l) s (2.37b) 

K(k) = P(k-l) c'^(k-l) (R+ C (k-1) P (k-l) C^(k-l))"^ (2.37c) 

P(k) = P(k-l) - K(k) C(k-l) P (k-1) (2.371) 

These equations are solved recursively given an initial value for the P 

A 

matrix, P(0), and an initial value for the estimates, q(0). 

21 

Mendel showed that if the expected value of the 
measurement noise is zero (i.e., E ( (k))= 0) and if and only if C(k) is 

free of measurement errors, then the parameter estimates will be unbiased in 
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the sense that 


E(q(k) ) = q 

The latter condition is clearly not applicable in the aircraft identification 
problem since C(k) is indeed constructed from noisy samples of the states. 
Then a biased estimate for q is euiticipated. 

Furthermore if the above condition did hold and if 
R, the vfei^ting matrix in 2.36 were selected to be the covariance of the 
actual measured value of 

x(k) - D(k-l) £ 

then q(k) as computed from 2.37 would be the linear minimum variance esti- 
mate for q, and P(k) would be the covariance for this estimate. 

Despite the inability to obtain a minimum variance 
using eq. 2.37# "tbe diagonal elements of R, r^^ , were in practice set equal 
to the variance in 

YiCk)-^ Di^ (k-1) s^ (2.38) 

0 


i.e. 





( 2 . 39 ) 


where j ranges over those indices for which D.,. ^ 0, and as in eq. 2.33 
2 

CT"^ is the variance of the disturbance inherent in y^ # the measurement of 
^i • 

To take into account the time variations in the 

2k 

parameters due to the aircraft maneuvers, either the fading factor cc can 
be incorporated or the equations can be rederived with the parameter variations 
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< 20 

modeled as fictitious noise sequences, ? ^ (k). Such a representation is: 

q (k+1) = q (k) + ^ (k), cov( f ) = Q (2.40) 

PracticaJLly speaking, this artifice causes the 
filter to continually track the parameters. The corresponding recursive 
scheme is similar to that of 2.37# hut now includes an update of the matrix . 


P(k) each sample, as given in: 

q(k) = q(k-l) - K(k) (C(k-l) q(k-l) - z(k) ) (2.1vla) 
z(k) = y(k) - D(k-l) s (2.4lh) 
P(k/k-l) = P(k-l) + Q (2.4lc) 
K(k) = P(k/k-l) C^(k-l) (R + C(k-l)P(k/k-l)C^(k-l) )"^ (2.4ld) 
P(k) = P(k/k-l) - K(k) C(k-l) P(k/k-l) (2.4le) 


The matrix Q has the effect of keeping P(k) from 
getting so small that parameter updating becomes insignificant due to a small 
gain matrix K(k). Rote that at steady state P(k) will go to a fixed value 
determined by Q. Selection of the elements of Q requires a trade-off 
between a filter that gives erroneous parameter estimates due to tracking of 
the measurement noise (if the elements of Q are too large) and a filter that 
just smooths out the parameter variations (if the elements of Q are too small). 

Computationally, the matrix inversion can be expected 
to cause trouble. However, if R and Q are chosen to be diagonal and if 
P(k) is initialized as a diagonal matrix, the term C(k-l)p(k/k-l)C (k-l) will 
be diagonal since C(k-l) contains only one nonzero- entry" in each colvimn. 

The corresponding matrix inversion is thus reduced to scalar divisions. 
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\ 

\ 


Details pertaining to the tuning of Q and P(0) 
are given in Section 4.2.1. 


2. 5 . 3 . 2 Stochastic Gradient Estimation 

In view of the fact that the weighted least squares 
algorithm ignores the noise present in the matrix C(k-l) in eq. (2.35) and 
consequently produces a biased estimate, a stochastic gradient procedure is 
being considered in anticipation of its capabilities for producing unbiased 
estimates. To illustrate its application, consider the i row of 


eq. 2 . 35 , i.e.. 




(2.42) 


To simplify the details it will be assimed that no teim in the i ^ row of 
C(k) is zero. Trivial modifications can be incorporated if some of these terms 
are in fact zero; 

Eq. 2.42 is now rewritten as: 


X. (k) - / D. .(k-l)s. 

1 Z_ 10 0 




(2.43) 


At this point it should be noted that since D. .(k-l) 
is either null, or is the state or control that couples x^ to s^ , a para- 
meter not requiring identification, the left hand side of (2.43) is a linear 
combination of measurable quantities. Consequently, a measurement z(k) for 
the left hand side of (2.43) can be constructed by replacing x^(k) and 
x^(k-l) by the respective corresponding noisy measurements y^|^(k) and 
y^(k-l). Thus 

z.(k) = X (k) - ^ D (k-l)s + (k) (2.44) 

X X . . XJ J 
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where Y" (k) contains the noise inherent in the measurement of x^(k) and 
in the measurement of those states necessary for foming I D (k-l)s . 


If C. .(k-l) in (2.1+3) is now replaced by its 
IQ 


corresponding measurement, then the following approximation model resxilts. 

X r..(k-l) q (2.1+5) 

ij J 


z,(k) 


where r. .(k-l) is the actual measured value for C..(k-1). 

Given this approximation model and assuming 
Tincorrelated noise sequences, the unbiased stochastic gradient algorithm for 


where: 


. 21 
is: 


^(k+1) = 

[l + W(k) ^ W(k)r^(k)Z^(k) 

^(k) is 

the estimate for q^. 



Z^(k) = : 

z-(k) - A r. .(k-l) q. . - . 1 . ' 

i' 10 0® 


( 2 . 1 + 6 ) 


r. (k) is the vector with components r. .(k) 
i 


I . 


is the covariance of the noise inherent to the 
measurement of r.(k-l) 


and . W(k) is a weighting matrix. 


To insure convergence W(k) must be of the form 
W(k) = ^ (k) diag |^h^(k), . . .,hj^ (k)"] 
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1+2 



Where: 


and 


0 -h ^h. - h < 0 D 

(. 1 a 

^(k) > 0 


OD 



k=o 


V ^ ^(k) < OD 

k=o 


These conditions essentially enable the procedure to he ahle to correct the 
parameter estimates as much as needed even for large values of k. 

Again because no simultaneous estimation of states 
and parameters is involved the algorithm is strictly linear. 


2,5.^ Procedure for Simultaneous State and Parameter Estimation 

Algorithms are now presented "ttiat perform both state and para- 
meter estimation. Those considered are the Extended Kalman filter and a 
decoupling process in \diich linear state estimation and linear parameter 
identification are performed separaitely and alternately. 

2 . 5.4.1 Extended Kalman Filter 

This section presents the general formula used in 
Kalman filter estimation without any attenpt for verification. Deteiiled 
derivations may be found in ref. (23) or any complete text dealing with 
filtering and estimation theory. 

The system defined by eqs. 2.32 and 2.33 vdll be 
used for the presentation of the identification algorithm. If only state 
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estimation were to be perfoimed, then the Kalman estimate for the state 

would be both a wii niTtnirn variance estimate and a least squares estimate in the 
sense of minimi zing the performance function^ 
n 

- H ip(k;)]’' [2p(R) - H^(k)] (2.47) 

k=0 

vrtiere R is the covariance of the measurement noise n. However if the 
measurement noise covariance matrix R is not known explicitly and/or liie 
noise is not Gaussian, then the estimate is optimal in the sense of minimizing 
n 

k=0 

when R is the estimated value of R. This estimate will not however be 
e 

optimal in the sense of minimizing the variance of the error. 

A , , 

The estimated state x^(k) is given by: 

^(k) = Ap Xp(k-l) + Bp np(k-l) + 

K(k) [yp(k) - H [Ap Xp(k-l) + Bp 

where: 

A 

Ap ®p "extrapolated" state. 

H l^p Xp(k-l) + Bp Up(k-lQ is the "extrapolated" observation, 

and K(k) is the gain matrix defined as: 

K(k) =[jd(k) [h M(k) + R^ (2.50) 

M(k) is the "extrapolated covariance matrix of 

X (k) given the observations y (k-l), y (k-2),... y (0) and is computed 
p P P P ■ 
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recursively as: 


M(k+1) = A P(k) A 


(2.51) 


P(k) is the covariance of the estimate x (k) given y (k), y (k-l),... 

P ' P P 

y^(0); that is 


r 


P(k) = E jp^Ck) - Xp(Er| [%(E) - 2Ep(k[|' 

and can be conrputed from M(k) and K(k) as: 

P(k) = M(k) - K(k) H M(k) 


yp(E)>...yp(o) 


(2.52) 


Up to this point, it has been demonstrated how Kalman filtering can be used 

A 

to estimate the plant state vector x^. 

To estimate any unknown parameters appearing 

in the state transition matrix Ap 'tbe gain matrix it is necessary 

to foim an augmented state vector. This is done by appending to the state 
equations the set of difference equations 

q^(k+l) = q^(k) i=l, ..., p (2.53) 

or in vector fpim 

q (k+1) = q (k) 

Thus the augmented system becomes: 

X®" (k+1) = A® (X^ (k), k) + B (X^k), k) u (k) 

y(k) = (k) X^ (k) + y\ (k) (2.54) 

where 

X = augmented state vector, given by 


X 


a- 


P 

q 
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.a 


A = matrix of the aupiented system given as 


a 

A = 






N, 


H = a\igmented system output matrix given hy 


H 


H 


N_ 


is an identity matrix of dimension ( p x p), and and are null 

matrices with dimensions of (n x p), (p x n) and (n x p) - respectively. 

Time variations in the parameters can be accoiinted 
for by including a fictitious noise component 2.53 such 


that: 


20 


q^(k+l) = q^(k) + ^^(k) 

where the variance and degree of correlation in are to be indicative 

of the anticipated fluctuation in q. A particularly useful assumption is to 
define ^ (k) as an •uncorrelated vector sequence with zero mean and covari- 
ance Q. This term would th--:n be added to the predicted covariance defined 
in eq. 2.51 giving: 

M(k+l) = A P (k) A*^ + Q 

An alternate method for tracking time varying para- 
meters is to introduce a fading on the filter's memory by replacing the 

25 26 

weighting matrix appearing in eq. 2.48 by: 


R(k) = a“^ Rg 


46 


a > 1 



where a the fading factor is typically between 1.0 and 1.1. 

Note that the augmented system, eq. 2.5U is non- 
linear and, therefore, incompatible, with the previous formula. This non- 
linear equation must therefore be linearized about a nominal or model 
trajectory x^. State vector and parameter estimates will then be given in 


terms of the incremental variations about this nominal trajectory or: 

\ 

K 

-P 

and 


X (k) = Xjj (k) + Ax (k) 


q (k) = qj^(k) + A q (k) 


(2.55) 

(2.56) 


where 3 ^ and q^^ are the state and parameters corresponding to the nominal 

or model system, and Ax (k) and A^ (k) are the corresponding incremental 

P 

quantities. Linearizing ( 2 . 5 ^) about the nominal augmented state vector 
T 

yields: 


A Xp (k+1) 
A q (k+1) 


J (k) 


where J (k) is the Jacobian matrix: 


J (k) = 


A (%(k)) 


A Xp (^) 
A q (k) 


(Ax + Bu) 




(2.57) 


(2,58) 


The extended Kalman filter equation applied to the 
( 23 ) 

system defined by (2.54) is thus given by' 
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AXp (k+l/k+1) 


t 

H 

A^ (k+l/k+1) 

= 

A q (k+l/k) 





vftiere: 


K(k) [AZp(k)-HAXp(k+l/k3 

( 2 . 59 ) 


A A 

AXp> are the estimated incremental variations in the state 

vector and parameters. 


Ay^(k) = (^) “ y^ 


(2.6o) 


and 


A Xp (k+l)/k 


Aq (k+l)/k 

Estimates for states and parameters at time k wotild thus he given by: 
^(k+l/fc) = x^(k+l) + A3^(k+l/k) 


= J (k) 


■ ^ / / s'* 

AXp (k/k) 


Aq (k/k) 


(2.6l) 


A 

q (k+i/k) 


% 


Aq(k+l/k) 


(2.62) 

(2.63) 


With regard to defining the model trajectory, the 
following possibilities should be considered: 

Procedure 1 ; 

The model or nominal trajectory can be defined by specifying the initial 
parameter vector qj^(O), and the initial state vector Xj^(O). These in turn 
define the solution to : 

y^^ (k) = H 3 ^ (k) (2.65) 
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where 


2^ = state vector of the nominal trajectory 
Ajj = system matrix defined hy the initial value of the 
parameters, ^(0). 

Bjj = gain vector defined hy the initial value of the parameters. 




, The matrix J(k) is then updated as follows using 
this nominal trajectory and the initial value of the parameters: 


J (k) = 



1 

I 

1 

1 

i 

1 

1 

% 

0 

1 

1 

1 

t 

1 



Procedure 2 t 

The nominal trajectory is defined using the previous estimate of the state to 
initialize the model equations at each step, i.e., 

Xj^(k+1) = Ajj Xp (k/k) + Bjj Up(k) 

As in the previous case 3^(k+l) is used to compute J(k) from eq. 2.58. 
Procedure 3 : 


The estimate q (k/k) is used to replace the model parameter vector q^, such 
that: 

Xj^(k+ 1 ) = A (q (k)) (k) + B (q (k)) u (k). 


and j(k) is subsequently computed using and q(k) as: 
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J(k) 


Ajj(q(k)) 


A 

0 

I 


Procedure 4 ; 

The equations are linearized about the most recent estimates for both 
and 3^. Thus 


3^(k+l) = Ap(q(k)) Xp (k/k) + Bp(q(k)) u^(k) 

(k+1) = q (k+l/k+1) 

This svimmarizes the necessary equations for the extended Kalman filter for 
State estimation and/or parameter identification. Details of the ntimerical 
experiments, condncted are presented in Section 4.2. 


2. 5. 4. 2 Separate Linear State and Parameter Estimation 
(SLSPE) 

The separate linear state and parameter estimator 

24 

( SLSPE ) as suggested by Graupe is a subop timal decoupling process in which 
linear state estimation and linear parameter identification are performed 
separately and alternately. Because each estimation process is basically a 
recursive wei^ted least squares algorithm, the procedures are basically those 
detailed in Section 2. 5. 3.1. 

To define the procedure consider again the linear system; 

X (k+l) = A X (k+l) + B u (k+l) (2.66a) 

-P P -P P -P 


Jp (1=) = 5p £ 


(2.66b) 
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when A and B contain xmloiown pEurameters. 

P P 

Define two linear systems, the state system defined hy eq. ( 2 . 66 a) and 
measured hy eq. 2.66h, and the parameter system: 

q (k+1) = q (k) (2.67) 


measured hy eq. 2.35- Clearly these two systems are coupled throu^ the 
system state - and the parameter vector 

In order to decouple these systems, states are to 
he estimated hy assigning some value to the parameter vector q and 
minimizing: 


N 


= I 


(y (k) - X (k) f E \ ^ 

P P * P P 


( 2 . 68 ) 


k=0 


subject to the satisfaction of eq. 2.66. 

Parameters are in turn to he estimated hy assigning a value to the state 


vector Xp and minimizing: 

N 


J = 5^ (y(J^) - x^(k) ) R (y(k) - x^(k) ) 

p *— • p p p 


(2.69) 


k=0 


where: 


Xp(k+1) = Ap(q) Xp(k) + Bp(q) u^(k) 


It is thus seen that the estimators for the states and 
the parameters are the algorithms defined respectively hy eq. 2.49 and eq. 

2 . 37 . In summary then, the SLSPE procedure is separable into the following 
two distinct problems: 
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A. At time k estimate x^(k/k), using for q 

its most recent estimate q(k-l^k-l) 

B. At time k estimate ^ (k/k), using for 

A 

its most recent estimate, i.e., x 

It is seen that the proposed procedure is suboptimal 
in the sense that the two objective functions are minimized separately and 
not simultaneously, resiolting in the need to alternate between a linear state 
estimation problem and a linear parameter estimation problem. 
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3.0 DIGITAL ADAPTIVE CONTROLLER REQUIREMENTS 


3.1 Considerations Concerning the Adaptive Algorithms 

When designing an algorithm to he implemented, practical considera- 
tions must he very influential on the trade-off between accuracy and sin^)!!-. 
city. In a digital adaptive flight controller, one of the prime practical 
restrictions is the size and speed of the on-line digital computer. This will 
affect hoth the timing and storage requirements of the adaptive algorithm. 
Figure 9 depicts the block diagram for an overall adaptive control system, and 
fig. 10 illustrates the various timing intervals involved in the realization. 

The smallest timing interval is the control computation interval 
(T^) and typically is between 0.03, and 0.20 seconds.. State estimation, if used, 
must be as fast, since the state is used by the feedback controller in calcu- 
lating the new control. The two larger intervals are required for gain 
adaptation. The parameter identification interval (T^) must not be longer 

than the gain update interval (T ), since the gain update algorithm requires 

S 

the new parameter estimates. If parameter estimation alone is perfoimed, 

then it is possible to have T^ > T^; however, if states and parameters are 

to be estimated simultaneously, then it is necessary that Gain 

update intervals (T ) being considered range from 1.0 to 10.0 seconds. 

S 

Since storage requirements vaiy from algorithm to algorithm, they 
must be taken into account. Recursive gain update, parameter identification 
and’ state estimation algorithms will require storage of the various matrices 


used in their formulation. 


Attention must be given to methods for reducing the complication of 
these algorithms. Prime among these is the identification of only the most 
sensitive parameters of the aircraft, the others being set at some representa- 
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FIG. 9 ADAPTIVE CONTROLLER STRUCTURE 
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Parameter 

identification 

interval 


FIG. 10 TIMING CHART FOR ADAPTIVE CONTROL 
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tive nominal value. The specific parameters to he identified are found 
according to the influence errors in identification have on the state. The 
sensitivity study performed for this puipose is described in Section 4.2. 

Furthermore, algorithms can be simplified by taking advantage of the 
sparseness of many of the matrices involved. This leads to a natural parti- 
tioning of the problem into a group or series of simpler problems. A further 
sin^jlification may be made for calculations involving symmetric matrices, 
although the extra information could be used to achieve more accurate calcula- 
tions . 

Structural considerations for implementing digital adaptive control 
systems based upon the single stage algorithm and linear quadratic regulator 
logic are next presented in Sections 3.2 and 3*3 respectively. 

3.2 Adaptive Logic Based Upon Single Stage Optimization Procedures 

Because of the characteristics outlined in Section 2.4.2, the con- 
troller based upon the single stage index vjhich penalizes the next stage error 
between plant and model was chosen for implementation in an adaptive control 
loop. 

The control law is derived based upon the minimization of a wei^ted 
combination of the instantaneous squared error between the states of a linear 
plant and model at the next stage and the present control energy. The result 
is a control signal that is a linear combination of plant and model states and 
model inputs as given in 2.4.2. The explicit foimula for computing the control 
gains are shown at the bottom of fig. 11. In order to satisfy the identifica- 
tion requirements, an on-line weighted least squares identifier, described by 
eq. 2.4l has been interfaced with this controller. The results of the 
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application of this controller to typical linearized lateral aircraft dynamics 
are discussed in 4.5.1. 

One of the prohlems encovmtered in the use of same state gains was 
that no set of gains gave acceptable response at any FC other than that for 
which it was designed, particularly for a sampling period of 0.2 sec. Thus 
in order to maintain stability and acceptable response with a same state con- 
troller it is necessary to update the controller gains regularly at a rate 
determined by the fastest esqpected rate of change of aircraft parameters in 
flight. 

With regard to implementation of the necessaiy logic for control 
computation, parameter estimation and gain update, the following operational 
procedures are required (assuming n states, m controls, 
meters to be identified in the i state equation ). 

CONTROL COMPUTATION 


and 


X 


para- 


Foimula: 


u (k) = K X (k) - K X (k) 

' X -in X — p' ' 

^ m P 


K 


u 


m 


u (k) 
-m' 


X (k+l) = A X (k) + B u (k) 

-in m -m' m -m 

Operational Requirements 

2 2 
Additions: 2m + 3nin - 3m - n + n 

2 2 

Multiplications; n + m + 3mn 


PARAMETER ESTIMATION 

Formula: eq. 2.4l 
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Opera tioneQ.. Requirements: ^ 

Additions; n(nH-n) + ^ 

i=l i=l 


n n 

Multiplications: n(nH-n+l) + ^ ^ '^i ^ 


i=l i=l 


n 

where ^ = total number of pareuaeters to be 


i=l 


identified. 


GAIN UPDATE 


Formula; 







II 

[h. 

b" 

p 

^®p- 


P 


II 

[- 

p 

Q®p. 




II 

[r. 

b’^ 

p 

^ V 





Operational Requirements: 
Additions : 


2 2 2 
3m + 2mn 3nm - 2m: + m 


2 2 2 

Multiplications; kan + mn - 2n 


Matrix inversions 
m = 1 
m = 2 


1 (m X m) matrix 
1 multiplication 

6 multiplications and 1 addition: 


As a specific example 
total computations would be necessary 


if n = 4, m = 2, and Jl^ = the following 
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Control: 


38 additions ^ 

> every T seconds 
44 multiplications^ ° 


Parameter identification: 

552 additions 

580 multiplication;^ 

Gain computation : 

83 additions 

134 multiplications j 


every seconds 
12 parameter estimated 


every T seconds 


For comparison, if an extended Kalman filter had been used for 
simultaneous state and parameter estimation, the above 552 additions would 
be replaced by 2235 additions and the 58 O multiplications would be replaced 
by 4580 multiplications. 

For T = T. = 0.2 sec. and T = 1.0 sec. the total additions and 

Cl g 

multiplications required per second would be: 

Additions: 3033 

Multiplications: 3254 

If only 4 parameters instead of 6 are identified in rowte 1 and 3, and 
= 0.1 second and = 1.0 second, the con 5 >utational requirements per 
sec. would be: 

Additions: 2043 

Multiplications: 2354 
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3-3 Adaptive Logic Based on the Linear Optimal Regulator Solution 

In view of the discussion in Section 2.4.3 an adaptive controller 

based on the solution to the linear optimal regulator problem was chosen for 

implementation as shown in fig. 12. Identification was performed, using the 

weighted least squares algorithm as defined in eq. 2.4l. Res\ilts using this 

adaptive logic are presented in Section 4.5.2. 

With regard to the implementation of the necessary logic for control 

con^nitation, parameter estimations, and gain update, the following operational 

procedures are required (assuming n states, m controls, andyf^ parameters 

til 

to be identified in the i state equation; 


CONTROL COMPUTATIONS 
Formula: 


Up(k+1) = Xp(k) + (I + K 22 ) u.p(k) + ujk) 

X (k+1) = A X (k) + B u (k) 
m' m mV m m' 


Operational Requirements: 

2 

Additions: 2m 

2 

Multiplications: 2m 


+ 3nin + 

+ 3™ + 


2 

n - m 



n 


PARAMETERS ESTIMATION 

■ Formula: Eq. 2.4l 

Operational Requirements: ^ 

Additions: n(m + n) + 2^^^^ + ^ 

i=l i=l 


n n 

Multiplications: n(m + n + l) + 

i=l i=l 
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PILOT 



ESTIMATES 


FIG.. 12 Digital Adaptive Controller Based bn Linear Optimal Control Theory 
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where 


n 


TJl. 


tot 


GAIN UPDATE 


Fonnula: (These were developed throiigh partitioning of the Riccati 

matrix, as defined in eq. 2.3I, talcing into account the 
sparseness of the A and B matrices as defined hy 
eq. 2.26). 


Pll 

P12 

P13 

Pl4 

T 

P12^ 

P22 

P23 

P24 

pis"^ 

P23*^ 

P33 

P34 

PlU"^ 

P24"^ 

p34'^ 

p44 


■pil(k+i) = Q + Ap'^(Pll - P12(P22 + R)‘^ P12'^) A^ (3.2a) 

P12(k+1) = Ap"^ [(PU - Pi2(P22 + R)"^ P12^) + P12(I - (P22 + R)"^ P22)] 

. • (3.2b) 

P22(k+1) = [(Pll - P12(P22 + R)'^ P12^) B^ + P12(I - (P22 + R)"^ P22)j 

+ [(I -• (P22 + R)'^ P22) P12^] Bp + P22(l - (P22 + R)"^^ P22) (3.2c) 

K22 = - (P22(k+1) + R)"^ (P12^(k+1) Bp + P22(k+1)) (3. 2d) 

= - (P22(k+1) + R)’^ (P12^(k-Hl) Ap) (3.2e) 

P13(k+1) = - Q + Ap'^(P13 - P12(P22 + R)"^ P23) A^ (3.3a) 
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P23(k+1) = Bp'^(P13 - P12(P22 + R)'^ P 23 ) A 

+ (P23 - P22(P22 + R)'^ P 23 ) A 

m 

^23 " " (P22(k+1) + R)"^ P23(k+1)A^ 


(3.3b) 


(3.3c) 


Pl4(k+1) - Ap'^[pi4 + (P13 . P12(P22 + R)'^ P23) - P12(P22+R)"^ P24 ] 

(3.3d) 

P24(k+1) = Bp’'[pi4 .* (P13 . P12(P22 ^ R)-l P23) 3^ . P12(P22^3J-1 pg^J 

(3.3e) 

+ (P23 - P22(P22 + R)"^ P23) B 

m 

+ P24 - P22(P22 + R)"^ P24 

^24 " " (P22(k+1) + R)-l (P23(k+1) + P24(k+1)) (3,2f) 


Operational Requirements; (per iteration) 

3 3 2 2 2 > 2 
Additions : 4n'^ + Sar + liftn n + m - I2mn - n - m 

3 3 2 2 

Multiplications: 4n'^ + 5m'^ + 9n m + l4m n 

Matrix Inversions: 1 (m x m) matrix 

m = 1 1 multiplication 

m = 2 6 multiplications and 1 addition 

As an example, if n = 4 and m = 2 the following operations are 


necessary: 


Control: 


42 additions 




48 multiplicationaJ 


every T seconds 


64 



Parameter identification: 

>f. = (4, 1, 3, 0) 

158 additions ^ 

178 raultiplications-J 


every seconds 


Gain computation: 

701 additions j 

*8l4 multiplications^ 


every T seconds 

o 


For = 0.1 sec. and = 1.0 sec. the operational requirements 

per second would he 

Additions: 27OI 

Multiplications: 3074 
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4.0 results 

4.1 Introduction 

Evaluation of a digital adaptive fli^t control system requires the 
testing of each individual subsystem prior to integration as well as escperi- 
ments on the entire system. Consequently considerable effort was expended in 
simulation efforts to individually evaluate each of the control algorithms 
discussed in Section 2.4 and each of the identification algorithms discussed 
in Section 2.5. 

Because of the need to examine the required preciseness of identifi- 
cation and the degree of axiaptation needed, a typical flight trajectory in the 
altitude-mach number plane was postulated. This is defined in Table II which 
cites the order and timing for a typical fighter aircraft to encounter ihe 
six given flight conditions of fig. 21. This trajectory corresponds to an 
initial acceleration from Mach .3 to ffech .9 at a very low altitude, a com- 
bined climb to 3600 Ja and acceleration to Mach 1.1, a climb to 15,000 m, a 
deceleration to Mach .9^ and finally a combined dive to 6000 m and a decelera- 
tion to Mach .7. For simulation purposes, it was assumed that the parameters 

of the aircraft's discrete A and B matrices varied linearly with time 

P P 

between flight conditions. 

A decision was made early in the study to use an adaptive control 
scheme based upon explicit identification of the unloiown parameters. However, 
as previously mentioned, it is impractical to consider estimation of all para- 
meters in the A and B matrices. A sensitivity study was therefore 
P P 

performed to determine those parameters which least effected system performance 
and which might therefore be considered constant. Results of this study are 
presented in Section 4.2, and resiilts of the identification experiments 
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Table II 


Flight Trajectory Used For Evaluation 


FC Time of encounter (seconds) 


1 

2 

3 

4 

5 

6 


0 

30 

35 

80 

85 

120 


6 ? 



themselves are presented in Section 4.3. 

Prior to testing the adaptive loop including the identifier, an 
evaluation of the gain update logic itself was made under the assumption of 
perfect identification of various influential sets of time vsuying parameters . 
This study, discussed in Section 4.4 led to conclusions regarding adaptation 
timing requirements and the number of parameters that had to be identified. 
Section 4.5 presents results for adaptive control incorporating online para- 
meter estimation. 

Finally, in order to assess the practicality of some of the control 
algorithms in a digital environment, several tests were made using NASA 
Langley's EAI 69 O hybrid computer. Both linear and nonlinear process dynamics 
were considered. These results are found in Section 4.6. 


4.2 State Sensitivity of Parameter Estimation Errors 
4.2.1 General Procedure 

The explicit calculation of the sensitivity vector 




S (k) 




(where q is some element of the discrete aircraft matrices) was used to rank 
the parameters as to the need for their identification. Depending upon the 
size of this sensitivity, it was decided either 

(1) to identify the parameter explicitly, thus having available an 
accurate estimate of its value, or 

(2) to set that parameter to some constant value and use that constant 

value for all gain update work. This constant (called was taken to be 

the average value of that parameter over all six flight conditions.. Thus, the 
error in the value of a parameter not Identified was Aq = q - 
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To get an approximate value for the error in. the state 
vector due to not identifying a given parameter, the parameter error was 
multiplied hy a representative value for the state sensitivity to that para- 


meter giving . 





z:^ 


<1 


( 4 . 2 ) 


• Thus even if. the sensitivity itself is large, a small parameter variation q 

over the flight envelope will yield a small ^ x^. Each parameter was then 

ranked, according to the error in each state of interest that would result from 

not identifying this parameter. . 

To calculate the sensitivity vector, a difference equation 

24 

was found that detennines the sensitivity vector's dynamics. This was done 
hy taking the partial derivative of the state equation: 


x(k+l) = f (x(k), u(k). A; B) = f(k) 
with respect to the specified parameter, q, resulting in: 


S (k.i) = lliii 

^ x(k) 


S(k) + 


^f(k) 


( 4 . 3 ) 


( 4 . 4 ) 


( 4 . 5 ) 


Since th-e sensitivity is of interest for the closed-loop system: 

X (k+1) = (A - B K) X (k) 

the sensitivity equations become: 

(A - B K) 

,S(k+l) = (A - B K) S(k) + X (k) 

■ ■ q 

Therefore, for a known state vector trajectory, the sensitivity of the various 
states to a certain parameter can be computed about that trajectory. It 
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shoTold be noted that these sensitivities will indeed be dependent upon both 
the system input signals u^ and the control gains implemented. 

Results of applying these sensitivity equations to the linear 
lateral equations of motion controlled by the single stage algorithms and the 

optimal regulator are next presented in Sections 4.2.2 and '4.2.3 respectively. 

• * 

4.2.2 Sensitivity Results for Single Stage Optimization 

Sensitivi'ty calculations for the controller based upon single 
stage optimization were performed using a simultaneous 5° aileron and 5° 
rudder pilot command input. In general, the set of most sensitive parameters 
for different flight conditions were the same except for some minor differ^ 
ences. The set of most sensitive parameters were obtained based upon the rank 
ordering of the error in roll rate and sideslip that would result if the 
average value of the parameters were used. It was found that in order to 
ensure good roll rate model following, the parameters that should be identi- 
fied, in order of importance are: B (1,2), A (1^4), B (l,l), A (l>2), 

p p B ■ P ' ■ 

Ap(1.3), Ap(l,l). 

To ensure good sideslip following the parameters in order of importance are: 
B_(3,2), A (3,4), B (3,1), A (3,2), A (3,3), A (3,1). 

x' ir jr , ir ir }r 

These twelve distinct parameters selected for identification consist of the 

elements of the 1st and 3rd rows of both A and B . Subsequent 

P P 

results reported in Sections 4.4.1 and 4.5.1 are based upon this set of 
parameters being identified. 
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4.2.3 Sensitivity Results for the Linear Optimal Regulator 

Sensitivity calculations for the opitmal linear regulator 
design were perfomed for a 5° aileron command acting alone and for a 5° 

; aileron command and 5° rudder command applied simultaneously. It was noted 
that, except for minor differences, the list of most important parameters was 
the same for all flight conditions. 

For the 5° aileron command, the following parameters 
(in order of descending in^jortance) should be identified to ensure good roll 
rate model following; 

Bp(l,l), Ap(l,l), Bp(l,2), Ap(3,l), Ap(l,2), 

To achieve accurate side slip response, the following parameters should be 
identified: 

Ap(3A)^ Bp(l,l), Bp(3,l), Ap(3^4) Bp(2,2). Note that there are ei^t 
parameters in the union of these two sets. 

It is also important to note that the list of most ir^ortant 
parameters is a function of the pilot input. For example, in comparison with 
the above and the results of Section 4.2.2, the 12 most sensitive parameters 
for a 5° aileron deflection and a 5° rudder deflection applied simultan- 
eously were: A (1,1), A (1,2), A (1,3)^ B_(l,l), B (1,2) 

.r }r h' 

Ajj( 2,3) Bp(2,2) 

Ap(3,l), Ap(3,2), Ap(3,3), Ap(3,4), B^(3>2) 

Note that this set differs from the single stage set of 4.2.2 in that it 
contains row 2 elements as well as some of the row 1 and row 3 elements . 
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4.3 Identification Results 

Evaluation of the various identification algorithms presented in 
Section 2.3 was perfonned through a series of experiments that enabled the 
appropriate design parameters to be tuned to the problem at hand prior to 
in^)lementing the fully adaptive control system. Considerations included: 
Performance in the presence of open loop control signals 
having relatively large bandwidths and variances. 

. Performance under open loop control signals having very 
little fluctuation. 

. Performance under closed loop control. 

Input requirements if elements of both columns of the 

matrix were to be estimated. For example if 

^ ^ = C ^ , then the corresponding control coupling 

to X . would be B (i,l)C + B (i,2)X = (B (i>l) +. 

B (i,2)C)^ ; and thus it would not be possible. to specify 

B (i.l) and B {x,2) uniquely. 

P P 

r Performance as a result of assuming erroneous measurement 
noise characteristics. 

Influence of the linearization procedures discussed in 
Section 2. 5.4.1. 

In general it was observed that the simple recursive wei^ted least 
squares algorithm gave better parameter estimates than did the extended 
Kalman filter. Apparently the improvement to be expected in simultaneous 
optimization of both parameter and state estimates is more than offset by the 
required linearizations. Consequently in subsequent adaptive control imple- 
mentations, the recursive weighted least squares algorithm was utilized. 
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Conments concerning the tuning of this algorithm are contained in Section 
4 . 3 . 1 , and examples of its perfoimance are contained in Section 4.5. As of 
this writing, conclusive results were not as yet available regarding 
the stochastic gradient algorithm, and stupes are being made of its 
usefulness in tracking time varying parameters in the presence of correlated 
noise. 

Perfoimance characteristics pertinent to the two algorithms for 
simultaneous state and parameter estimation, i.e., the extended Kalman filter 
and the separate linear state and parameter estimation algorithms (SLSPE) are 
contained in Section 4.3*2. 

4 . 3 .I Results Pertinent to Procedures for Parameter Estimation 

Only 

The study conducted thus far on the stochastic gradient 
estimation has been basically a theoretical one and has been geared towards 
the development of the procedure for the identification of constant parameter 
multi-input multi-output systems. Further development is expected in the near 
future to adapt the procedure to the identification of time varying parameters 
and for the establishment of experimental results and analysis. This research 
will be concentrated on the synthesis of the correction gain W(k) given in 
eq. 2.46, taking into account both simplicity and speed of convergence, and on 
the 'evaluation of performance in the presence of correlated noise sequences. 

With regard to the recursive weighted least squares algorithm, 
some of the results pertinent to the selection of the P(0) and Q matrices 
as used in eq. 2.4l are now presented. As stated in Section 4.1, 
this algorithm was selected for implementation in the fully adaptive system. 
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The pertinent parameter tracking capabilities are thus discussed in Section 
4.5 eilong with the other results obtained from testing the overall system. 

The use of large values in the Q matrix results in the 
ability to track rapidly varying parameters. However, it also results in a 
tendency to adjust parameters so as to account for variations due solely to 
measurement noise. Similarly, small values in the Q matrix result in a more 
pronounced filtering of the measurement noise while the ability to track time 
varying parameters is degraded. 

As a first approximation, the appropriate elements of Q were 
set equal to the expected square of the appropriate parameter change per 
sample, and P(0) was then set equal to the corresponding steady state solutions 
of eqs. 2.4lc, 2.4ld, and 2.4le. In practice, this was accomplished by iter- 
ating throu^ these filter equations and observing the steady state value of 

P. Elements of the R matrix were computed in accordance with eq. 2.39, 

• 2 

although it was observed iiiat setting r^^ equal to gave equally good 

resTilts. 

For estimation of the 12 most sensitivity parameters as used 
in the single stage adaptive algorithm (see Section 4.2.2), the diagonal 
elements (all other elements were zero) of P(0) and Q were: 
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Parameter 


P. .(0) 
ix' ' 





8.0 X 10"^ 

4.?r X 10“^ 

yi,2) 

2.0 X 10“3 

1.1 X 10”^ 

yi,3)' 

1.6 X 10"^ 

7.75 X lO"^ 

yi,4) 

6.5 X lO"^ 

6.4 X 10 

ys,!) 

1.4 X 10“^ 

4.4 X 10“^ 

Ap(3,2) 

2.9 X 10'^ 

4.0 X 10"® 

Ap(3,3) 

2.5 X 10“^ 

7.1 X 10"® 

V3,4) 

1.4 X 10"^ 

4.4 X 10“^ 

Bp(l,l) 

1.9 X lO"^ 

9.4 X 10"^ 

Bp(l,2) 

4.4 X 10“^ 

6.4 X lO"^ 

Bp(3,l) 

3.8 X 10’^ 

1.0 X 10"® 

Bp(3,2) 

5.9 X 10“^ 

1.1 X lO""^ 


For the eight most sensitive parameters as used in the adap- 
tive controller based upon optimal regulator logic (see Section 4.2.3) the 


elements of P( 0) 

'and Q were: 


Parameter 

P. .(0) 
ii 



2.0 X 10"3 

6.0 X lo"^ 

Ap(l,2) 

2.0 X 10"3 

1.6 X lO""^ 

B^(IA) 

1.2 X 10"^ 

3.4 X 10"^ 

Bp(l,2) 

5.0 X 10"^ 

-4 

2.3 X 10 ^ 

Bp(2,2) 

5.0 X 10"^ 

3.4 X 10"® 

y3,i) 

1.2 X 10"^ 

. 8.6 X 10"^® 

V3,4) 

1.5 X 10'^ 

1.2 X 10"'^ 

Bp(3,l) 

-4 

1.5 X 10 ^ 

c 

2.0 X 10 
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4.3.2 Results pertinent to Estimation of Both States and Parameters 
Conclusions pertinent to the perfo rman ce characteristics of 
both the Kalman filter and the SLSPE were based upon simulation on an IBM 
360/50 digital con^juter. Sampling periods of I/30 second and l/lO second were 
incorporated, and estimation of the four lateral states and the eight most 
sensitive parameters for the linear regulator logic. (as presented in Section 
4.2.3) was considered. Unless otherwise cited, measurement noise was simulated 
as uncorrelated processes having the variances given in Table I. 

4. 3.2.1 Kalman Filter RestLLts 

The following conclusions were reached regarding the 
perfoimance of the extended Kalman filter: 

Selection of initial covariance matrix P(0) 

Assuming perfect knowledge of the initial state, the 
corresponding elements of P(0) should be zero. 

The terms of P(0) corresponding to the parameters 
can be zero if fictitious noise with a covariance 
, matrix Q is used to model the parameter variation; 

otherwise these elements should be proportional to 
the expected square of the per sample change in the 
parameter. 

Measurement noise covariance matrix 

It was noted that the use of the exact noise covariance 
in the filter equations was not necessary for con- 
vergence. It was observed that convergence was possible 
with the noise covariance R set five times larger than 
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the actual noise covariance matrix; however 
convergence did not occur for R matrix element-S' 
much lower than those in the actual measurement 
noise covariance matrix. 

Nominal trajectory selection 

For time varying parameters it was foiuid that the 
hest' procedure for updating the nominal trajectory 
was throu^ the use of the most recent state and 
parameter estimates for defining the point about 
which linearization was to be performed. This is 
Procedure 4 as discussed in Section 2. 5.4.1. For 
constant parameters, it was found that the best 
trajectory to linearize about .was defined by 
Procedure 2. 

Fading memory implementation 

From experimental results, it was observed in general 
that the use of a fictitious noise input for modeling 
the parameter variations yielded better convergence 
properties than did the procedure using an experi- 
mental weighting factor. This is illustrated in 
fig. 13 for estimation of a^^^ under open-loop 
excitation (+ 30° rudder at I .5 Hz and + 30° aileron 
at 1.0 Hz). State measurement noise was uncorrelated 
with variances of one tenth those listed in Table I. 
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FIG. 13 a TYPICAL PARAMETER TRACKING USING FICTITIOUS NOISE 
TO MODEL PARAMETER VARIATIONS 


• • • ESTIMATES 
TRUE VALUES 



FIG. 13 b TYPICAL PARAMETER TRACKING USING A FADING FACTOR 
IN THE PERFORMANCE INDEX 
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4. 3. 2. 2 SLSPE Results 


In using the SLSPE it was observed that the para- 
meters estimates were biased and that convergence was poor if parameters in 
both columns of the matrix were identified. However, these results 

corresponded to experiments in which a recursive step of both the parameter 
estimation algorithm and the state estimation algorithm were performed at each 
and every sairple. In view of the restilts obtained using the weighted least 
squares algorithm for parameter estimation only, it is recommended that 
either state estimation not be performed as often as the parameter estimation 
or that the state measurements, not the state estimate^ be used in the para- 
meter estimation algorithm. Future efforts w ill consider these items. 


4.4 Performance Under Perfect Identification Prior to Interfacing 

Adaptive Control Logic 

With an online identifier, it is necessary to validate the results 
of the sensitivity study described in 4.2.1 and to deteiming how often adapta- 
tion shotild be performed. 'Upper bovinds to these figures can be obtained by 
simulating an adaptive controller with a perfect identifier. Thus while 
parameters not being "identified" would be assumed at their average values, 
true values would be used for these parameters . designated for identification. 

4.4.1 Performance of the Single Stage Algorithm 

The following two experiments pertain to the performance of 
the' single stage controller: 

Experiment I 

Objective: Assuming no noise and perfect identification of all para- 

meters, determine suitable gain update periods. 
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Procedure: At each gain update sample time, the gains optimal for the 

current flight conditions were to be computed from eq. 2.19 and held constant 
until the next gain update sample time. These gains at the 6 PC's are given 
in Appendix B. 

Constant factors: . 

Pilot input: + 5° aileron, 0.1 hz square wave 

: Variable factors: 

Control sampi i ng periods: .2, .033 sec. 

Gain update periods: 1.0,, 2.5, 5*0 sec. 

Results: Figure l4 indicates good model following for a gain update 

period of 2.5 sec. and a control period of 0.2 sec. A gain update period of 
5.0 sec. was found to be satisfactory for a control san^jling interval of 
.033 sec., but \ansatisfactory for a control period of 0.2 sec. 

Experiment II 

Objective: Assuming that only selected sets of 10-12 parameters are to 

be identified perfectly, with the remainder set equal to their avereige value 
over the typical flight trajectory for six FC's determine the necessary 
frequency of gain updating. The parameter sets to be identified are those 
deemed to be the most influential as discussed in 4.2.1. 

Procedure: At each gain update sample, the gains were to be computed from 

eq. 2.19 using the true values for those. parameters being "identified" and 
average (over the 6, FC's) values for these parameters not being identified. 

Constat factors: 

Pilot input: + 5° aileron, 0.1 hz square wave 

Control sample period: 0.2 sec.. 
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Variable factors: 


Parameters identified: 10 and 12 most influential parameters 

Gain update periods: 1.0, 2.5 sec. 

Results: Figures 15 and l6 indicate the significant improvement that 

results from identifying 12 rather than 10 parameters especially between 
FC3 and FC4. This is true for both gain update periods. The most notable 
improvement in using a gain update period of 1.0 rather than 2.5 sec. is most 
evident between FCl and FC2. 

Conclusions of both e3q)eriments : For the overall adaptive implementation, 

identify the 12 most influential parameters, and based upon these, update 
gains every 1.0 sec. 

4.4.2 Performance of the Linear Optimal Regulator 

The following three experiments pertain to the perfomance 
of the linear optimal regulator design: 

Experiment I 

Objective: Assuming no noise and perfect identification of all para- 

meters, determine s\iitable gain update periods. 

Procedure: At each gain update sample time, the P matrix is recomputed 

from a single iteration of the Riccati eq. (2.31) using those values of A^ 
and Bp which correspond to the present flight conditions. Control gedns 
are then computed from eq. 2.29. (See Appendix B for the gains corresponding 
to the 6 given PC's). 

Constant factors: 

Pilot input: + 5° aileron, .2 hz square 
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FIG. 15 SINGLE STAGE ADAPTIVE CONTROL LOGIC- ROLL RATE RESPONSE 
CONTROL INTERVAL 0.2 SEC 
GAIN UPDATE INTERVAL 1.0 SEC 

PERFECT IDENTIFICATION OF 10 AND 12 PARAMETERS 
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ROLL RATE (DEG/SEC) 



0 2 4 6 8 10 30 32 34 36 74 76 78 80 82 84 


TIME (SEC) 

FIG. 16 SINGLE STAGE ADAPTIVE CONTROL LOGIC -ROLL RATE RESPONSE 
CONTROL INTERVAL 0.2 SEC 
GAIN UPDATE INTERVAL 2.5 SEC 

PERFECT IDENTIFICATION OF 10 AND 12 PARAMETERS 


84 


Control sample period: 0.1 sec. 

Variable factors: 

Gain update periods: 1, 5, 10 sec. 

Results: Figures 17a and b show good model following for a gain update 

period of 1.0 sec., while figs. l8a and b indicate some degradation in perform- 
ance for a gain update period of 5.0 sec. Stability for the closed loop 
system was observed for the three tested gain update periods of 1, 5 and 
10 sec. 

Experiment II 

Objective: Under the assmption of no noise and perfect identification of 

the 8 parameters, chosen by the sensitivity method described in Section 4.2, 
with all others set to their average value over all 6 flight conditions, 
determine the required frequency of gain update. 

Procedure: At each gain update sample, the gains were to be computed 

from one iteration of the Riccati eq. (2.31) using the true values for the 
ei^t parameters being identified and average (over the 6 FC's) values for 
the remaining l6 parameters . 

Constant factors: 

Pilot input: + 5° aileron, .2 hz square wave 

Control sample period: 0.1 sec. 

Variable factors: 

Gain update periods: 1, 5^ 10 sec. 

Results: Figs. 19 and 20 illustrate the improvement in model following which 
results from using a 1 sec. gain update period rather than a 5 sec. period. 

For a 10 sec. period, the stability of the closed loop system was marginal. 
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FIG 17a ROLL RATE RESPONSE 



5 10 30 35 40 75 80 85 

TIME (SEC) 


CONTROL INTERVAL 0.1 SEC ' 

GAIN UPDATE INTERVAL 1.0 SEC — MODEL ' 

PERFECT IDENTIFICATION OF ALL PARAMETERS — AIRCRAFT 



5 10 30 35 40 75 80 85 

TIME (SEC) 


FIG. 17 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC- 
ROLL RATE AND SIDESLIP RESPONSES TO 
i5®, 0.2 hz AILERON COMMAND 
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FIG. 18a . ROLL RATE RESPONSE 



O 5 10 30 35 4 0 75 8 0 85 


TIME (SEC.) 


CONTROL INTERVAL 0 1 SEC 
GAIN UPDATE INTERVAL 5 0 SEC. 
PERFECT IDENTIFICATION OF 
ALL PARAMETERS 


MODEL 

AIRCRAFT 



0 5 10 30 35 4 0 75 80 85 

TIME (SEC) 


FIG. 18 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC- 
ROLL RATE AND SIDESLIP RESPONSES TO 
0.2 hz AILERON COI\JlMAND 



FIG. 19 a ROLL RATE RESPONSE 



O 5 10 30 35 40 75 80 85 


TIME (SEC) 


CONTROL INTERVAL 0.1 SEC MODEL 

GAIN UPDATE INTERVAL 1.0 SEC AIRCRAFT 

PERFECT IDENTIFICATION OF EIGHT PARAMETERS 



TIME (SEC) 


FIG. 19 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC - 
ROLL RATE AND SIDESLIP RESPONSES TO 
^5°, 0.2 hz AILERON COMMAND 
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I FIG. 20a ROLL RATE RESPONSE 



O 5 10 30 35 40 75 80 85 


TIME (SEC) 


CONTROL INTERVAL 0,1 SEC 

GAIN UPDATE INTERVAL 5.0 SEC 

PERFECT IDENTIFICATION OF EIGHT PARAMETERS 


MODEL 

AIRCRAFT 



FIG. 20 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC - 
ROLL RATE AND SIDESLIP RESPONSES TO 
^5°, 0.2 hz AILERON COMMAND 
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Note that relative to the results of Experiment I, there is more deterioration 
of side slip model following than there is of roll rate following. This is 
true because the set of eight parameters identified is more influential in 
defining the roll rate dynamics. 

Experiment m 

Objective: Deteimine how the results of Experiment II are affected if 

noisy state measurements are used for control computation. 

Procedure: Same as for Experiment II except control signals were to be 

computed from products of the gains and noisy state measurements. 

Constant factors : 

Pilot input: + 5° aileron, .2 hz square wave 

Control sample period: 0.1 sec. 

Gain update periods: 1, 5 secs. 

Variable factors: 

Noise statistics: Table I characteristics achieved using 

first order filters with breaks at the 
given bandwidths. RMS values given in 
Table I and l/lO these values were used. 
Results: Figure 21 shows that feeding back states contaminated with 

noise results in severe Jittering of the aircraft response. Some sort of 
prefiltering is thus necessary. As expected, lowering the rms values of the 
noise samples reduced the magnitude of this Jitter. 

Conclusions from all three experiments: For implementation of an overall 

adaptive system, identify the eight most influential parameters discussed in 
Section 4.2.3, andusing these, update the gains every one to five seconds. 
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0 5 10 30 35 40 75 80 85 

TIME (SEC.) 


CONTROL INTERVAL 0,1 SEC 
GAIN UPDATE INTERVAL 5 0 SEC. 

PERFECT IDENTIFICATION OF EIGHT PARAMETERS 
FEEDBACK OF ACTUAL MEASUREMENTS 


MODEL 

AIRCRAFT 



0 5 10 30 35 40 75 80 


TIME (SEC.) 

FIG. 21 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC - 
ROLL RATE AND SIDESLIP RESPONSES TO 
i5°, 0,2 hz AILERON COMMAND 
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Prefiltering is needed for state measurements contaminated 'by noise having 
the characteristics of Table I. 

4.5 Adaptive Control System Performance Evaluation 

4 . 5.1 Adaptive Control Loop E 3 q)eriment - Single Stage 
Optimization Procedure 

Objective: To study the behavior of the adaptive control system designed 

by coupling a weighted least squares identifier with the single stage con- 
troller logic as shown in fig. 11. 

Procedure: Using noisy state measurements with no prefiltering, para- 

meter estimates were obtained at each control sample period and then used at 
each gain update period to compute the control gains according to eq. 2.19. 

The noisy state measurements were then used with these gains for control 
conputation . 

Design factors: 

Pilot input: + 5 ° aileron, 0.1 hz square wave. 

Control sanple period: 0.2 sec. . 

Gain adaptation period: 1.0 sec. 

Parameters identified: Ap(l>l)^ Ap(l,2), Ap(l,3)^ Ap(l^^)# 

Ap(3,l)> Ap(3,2), Ap(3,3), A^(3,4), 

Bp(l,l), Bp(l,2), Bp(3,l), Bp(3,2). 

Remaining parameters: Set at the average values' computed over the 6 given 

PC's. 

Measurement noise: As defined in Table I. 
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Resiolts: 

Responses 

Figure 22 shows the roll rate response for portions of 
the test flight trajectory. The noisy state feedback results 
in severe jittering of the aircraft and thus relatively poor 
model following particularly at FC's 1 and 2. Side slip was 
generally limited to + 0.4° except for FCl where the maximum 

t 

excursion was 2.8°. 

Parameter tracking 

Figures 23a, h, c, 'd i3JLustrate parameter tracking 
"behavior between FCl and FC2 during which model following is 
the worst. To he noted is that despite the existence of a 
bias in the estimate, the general trend in the variation is 
followed. 

Control gains 

Comparisons between the optimal and adapted control gains 
are given in figs. 24a, b, c, d for part of the interval 
between FCl and FC2. With the exception of the roll rate 
feedback gains, the adapted gains were relatively close to 
the optimal gains. In fact, beyond the 20 second interval 
depicted in fig. 24, all gains converged closer to the true 
• values. 

Conclusions: In view of the fact that the adaptation procedure was capable 

of producing reasonable feedback gains, it is concluded that a workable system 
should res-ult from the inclusion of prefilters for removing jitter and the 
usage of an xmbiased parameter identification algorithm. 
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FIG. 22 ROLL RATE RESPONSE, SINGLE STAGE LOGIC 

WEIGHTED LEAST SQUARES IDENTIFICATION OF TWELVE PARAMETERS 
CONTROL INTERVAL 0.2 SEC 
GAIN UPDATE 1.0 SEC 
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4.5.2 Adaptive Control Loop Experiment - Linear . Optimal 
Regulator Procedure 

Objective: To study the behavior of the adaptive control system designed 

by coupling a weighted least squares identifier with the Riccati gain update 
logic as shown in fig . 12 . 

p2X)cedure: Using noisy state measurements with no prefiltering, para- 

meter estimates were obtained at each control sample period and then used at 
each gain update sample in one iteration of the Riccati eq. (2.31). The gains 
were then computed from eq. 2.29, and along with the noisy state measurements 
used for control computation. 

Design factors: 

Pilot input: + 5° aileron, .2 hz square wave 

Control san^le period: 0.1 sec. 

Gain adaptation periods; 1, 5 sec. 

Parameters identified: Ap(l,l), A^(l,2), A^(3,l), Ap(3,4), 

,B (1,1), B (1,2), B (2,2), B (3,1). : 

Jr Jr x' 

Remaining parameter;. Set at the average value computed over the 6 given 
PC's.. 

Measurement noise: As defined in Table I. 

Results : 

Responses 

Figure 25a and b depict roll rate and side slip responses 
for a gain update of ’^.0 sec. These responses are very 
similar to those of fig. 21 (for noisy feedback but with 
perfect identification of the above 8 parameters) indicating , 
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0 5 10 30 35 40 75 80 85 

TIME (SEC.) 


CONTROL INTERVAL 0.1 SEC. 
gain UPDATE 5.0 SEC. 

WEIGHTED LEAST SQUARES IDENTIFICATION OF 
EIGHT PARAMETERS 


MODEL 

aircraft 


FIG. 25b SIDESLIP RESPONSE 



TIME (SEC.) 

FIG. 25 LINEAR REGULATOR ADAPTIVE CONTROL LOGIC- 
ROLL RATE AND SIDESLIP RESPONSES TO 
i5®, 0.2 bz AILERON COMMAND 



that the inclusion of an actual identifier did not seriously 
degrade the performance. With a gain update period of 1.0 
sec., improvements in the model following chaxacteristics 
were noted. 

Parameter tracking 

Figures 26a, h and c illustrate typical parameter 
tracking. Again as in ^+.5.1 the presence of bias in the 
estimate is readily evident as well as the capability for 
tracking the variation. As might be escpected, a reduction 
of the nns values to l/lO those cited in Table I improved 
the estimates. 

Control gains 

Con^arisons between the gains computed based upon 
perfect identification of all parameters and the actual 
adapted gains are shown in figs. 27 a, b, c, d, for part of 
the interval between FCl and FCU, Despite the differences 
between the adapted and optimal gains (most notably in liie 
feedback gains affecting roll rate) the general variations 
in the gains were followed. 

Conclusions : In view of the facts that the procedure was capable of adapting 

the feedback gains to reasonable values and the simularity between figs . 21 
and 25 , it is expected that the incorporation of prefiltering and an unbiased 
parameter estimation algorithm should produce an acceptable adaptive control 
system. 
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FB FROM XP(3)-SIDESLIP FB FROM XP (l)-ROLL RATE 



FIG, 27 FEEDBACK. CONTROL GAINS - LINEAR 


FB FROM XP(4)-BANK ANCLE ^ pg pp,Qj^ XP(2)-YAW RATE 
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4.6 I^brid Cpioputer Tests 

In addition to the previously described 1 digital simulations^ 
the EAI 690 analog- digital hybrid conputer at NASA LRC was used for system 
evaluation. The aircraft dynam i cs and pilot controls were programmed on the 
analog portion of the computer, while the discrete control laws were pro,- 
grammed on the digital portion. 

4.6.1 Linear Simulation Results 

The first tests were perfoimed using the linearized aircraft., 
dynamics. These tests had as their objectives: 

Evaluation of feedback controllers in a hybrid environment. 

Evaluation of the gain adaptation algorithm using perfect 
parameter identification. 

The restilts indicated that: 

Noise in the analog to digital converters can be a significant 
problem for the nondynamic feedback controller developed using the 
single sta^e algorithm. The "type one" feedback controller developed 
using linear optimal control theory was much less sensitive to this 
noise. In particular, a noisy side slip measurement caused trouble 
for the single stage controller because of the moderately large 
gains necessary to effect good model following. 

Even when no information about the actual flight condition was 
included in the initial Riccati matrices, only one or two executions 
of the gain update algorithm were necessary to produce gains giving 
accurate model-following. Figure 28a illustrates the aircraft 
response to a pilot aileron input using incorrect gains (those 
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designed for FCl actually applied to the aircraft at FCS), while 
fig. 2Sb illustrates the aircraft response after just one iteration 
of the gain update algorithm. This aircraft response is indistin- 
guishable from the model response. 

4.6.2 Nonlinear Simulation Results 

A second set of preliminary tests were performed using a 
nonlinear 6 degree-of-freedom aircraft analog simulation and the same d^amic 
feedback controller as described above. To avoid the logic needed to account 
for discontinuities in aircraft bank angle at 0 = l 80 °, the sine, of this angle 
was fed back. 

The lateral control system maintained the same performance 
for small perturbations about the trim condition as it had in the linear case, 
and did not fail to operate under extreme maneuvers . Figure 29 illustrates 
some of this behavior. 
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SIMULATION OF A FULL ROLL 



5.0 DISCaSSION MD CONCLUSIONS 

5.1 Results rdscussion 

In summary, the results presented in Section 4.0 lead to the follow- 
ing observations : 

. Digital adaptive flight controllers can be designed using either 
linear- quadratic optimal regulator theory or single stage perfonn- 
ance functions . 

. Because of superior stability margins, the design based upon optimal 
regulator theory requires fewer parameters for identification. 

. The perfoiTnance of the overall adaptive control system fig. 25 is 
comparable to the perfoimance of the "adaptive" controller with 
perfect identification simulated fig. 21. This was based upon 
ejcperiments using measurement noise with the characteristics of 
Table I . 

. The adapted control gains in general track the variations in the 
optimal gains. However, some bias is evident especially in the 
feedback gains involving roll rate. 

. Prefiltering of the noisy measurements prior to use in feedback 
control computation is needed. 

. Usage of the noisy measurements themselves in the weighted least 
squares identification algorithms appears permissible. 

. Identification of only parameters, as opposed to estimation of 
states and parameter together, appears preferable in view of the 
results obtained and the computational and linearization require- 
ments of the latter. 
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Application of the linearly designed controller to the actual 
six degree of freedom coupled system should yield acceptable 
' - results in view of the results discussed in Section 4.6.2. 

5.2 Feasibility of Implementation 

Both of the adaptive control algorithms discussed in the preceding 
should be re&dily in: 5 )lementable in a digital flight computer. For purposes 
of illustration, consider an add time of three micro-seconds and a multiply 
time of six micro -seconds. 

For the single stage design discussed in Section 3-2 the following 
timing requirements would be typical for a fourth order system with two 
control inputs and with weighted least squares identification being performed 
for 12 parameters; 

Control computation .378 ms 

Parameter identification 5.13 ms 

Gain confutation I.05 ms 

For estimation using an extended Kalman filter the timing requirement would be 
34.2 ms. Thus if at the worst the above three functions were performed at 
each and every sample, a minimum sample period of about 7 ms would be needed 
for weighted least squares identification and about 4o ms if extended Kalman 
filtering were used. 

Similarly for the optimal design described in Section 3*3^ the 
following woiild be typical timing requirements for a fourth order system with 
two controls and weighted least squares • estimation of eight parameters 
(identification of 12 parameters was not necessary): 

. 4l4 ms 
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Control confutation 



Parameter identification 


1 . 54 ms 

Gain con5)utation 6.99 ms 

Thus if in the worst case these three computations had to be performed at 
each and every sample, a minimum period of about 10 ms would be required. 

5.3 Recoimnendations for future research 

In view of the experimental results discussed in Section 4.0, it is 
recommended that the following topics be studied: 

. Design of filters to be used for smoothing the measurements 
prior to control computation. 

The use of prefiltered measurements for identification. 

. A stability analysis of the single stage design perhaps 

. utilizing the procedures followed by Chan^^. These would • 
incorporate the use of a perfect model following controller 
despite the non-satisfaction of the conditions of perfect . 
model following . 

. Applicability of the linear controller design to the o 
nonlinear six-degree-of freedom system which incorporates 
coupling between lateral and longitudinal motion. 

Performance when true actuator dynamics, bending modes, 
and turbulence are considered. 

. The need. for dither for improving identification. 

. The relative performance of an implicit adaptive controller 

based on hyperstability principles (see e.g., Landau, ref. 6). 
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Appendix A 


Aircraft Continuous LateraJ. Matrices 
Corresponding to Equation 1.1. 
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Appendix B 


SS Lateral Controller Gain Matrices, 

Q = diag (1, 0, 1, 0), R = diag (0, 0), Sampling 
Period .2 sec. 
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Appendix B 


Optimal Regulator Gain Matrices, 

Q = diag (lo\ 0, lo\ 0), R = diag (10^, 10^), 
Sampling Period .1 sec. 


FCl 



P r 

P 


II 

H 

<M 

-.302 - 1.04 

4.4 

.081 

-.145 1.42 

-1.7 - 

.II5J 


Feedback 

from X 




P 


K 23 

f.0975 .55 

- 1.10 

0.0I 

1_.025 -1.15 

1.25 

o.oj 


Feedback 

from X 
m 


FC2 




II 

-.025 -.50 

1.7 

.021 

-.005 +.47 

0.3 - 

.02J 

K 23 - 

r.025 ,35 

- 1.0 

0.0] 

L.OO75 -.70 

1.75 

o.qj 

FC3 




^1 “ 

-.040 -.67 

2.0 

.025I 

-.010 .45 

0.8 - 

.015J 

K 23 = 

f.0275 .53 

-1.3 

0.0] 

[.010 -.75 

1.65 

O.OJ 



- 1.54 - .08] 

- .30 -I.28J Sj ' 


Feedback from u 
. P 



Feedback from u 


K22 = 
K24 = 

K22 = 
K24 = 



114 



fc 4 


K21 = 

-.54 

L -23 

- 1.2 

1.6 

4.9 

- .6 

.035^ 

-.05 

K 22 = 

F1.63 

L.37 

-.20 I 
-I.06J 

K23 = 

rn 3 

L0275 

.58 

-1.25 

-•95 

0.90 

o.o" 

0 . 0 _ 

k ^4 = 

fl .53 

L -13 

-.281 

I.20J 

FC5 ' 

p 

r 

P 





^ 21 = [ 

-.54 

r.23 

-1.25 

1.80 

5.1 

-1.3 

.0451 

-.065J 

K22 = 

-1.66 
r -37 

-.23 1 
-I.I8J 

K 23 = 

~.108 

_.035 

.63 

-1.25 

-.95 

• 6> 

o.o" 

0 . 0 _ 

% = 

'1.52 
- -15 

-.34 1 
1 . 28 J 

FC6 








K21 = 

-.18 

-.07 

-.85 

1.15 

3.7 

-1.3 

.o45~ 

-. 055 . 

K22 = 

I1.56 

r -23 

o.o” 

-i.4_ 

K23 = 

r.058 

L.oi5 

.45 

-1.00 

-1.0 

1.5 

0.0' 

0 . 0 _ 

K 24 = 

r.75 

L-13 

-. 22 " 

.74_ 


'itv. S. GOVERNMENT PRINTING OFFICE: 1974-640046/44 


115 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE S300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND PEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
481 



POSTMASTER : 


If Undeliverable (Section 168 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere artd space. The Administration 
shall ^ovide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof." 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES; Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA aaivities. 
'Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTIUZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORAAATION OFFICE 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



